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Abstract 



The long life-time of the quark-gluon plasma likely to be created in the relativistic 
heavy ion collisions at RHIC and LHC energies renders it sensitive to the details of 
the equation of state of the quark-matter. We show that the pr dependence of the 
survival probability of the directly produced J/ip at RHIC energies and that of the 
directly produced T at LHC energies is quite sensitive to the speed of sound in the 
quark matter, which relates the pressure and the energy density of the plasma. The 
transverse expansion of the plasma is shown to strongly affect the J/tp suppression 
at LHC energies. 



PACS: 12.38.M 



1 Introduction 



Statistical quantum chromodynamics predicts that at sufficiently high densities or 
temperatures the quarks and gluons confined inside hadrons undergo a deconfining 
phase transition to a plasma of quarks and gluons. The last two decades of high 
energy nuclear physics activity has been directed towards the production of this new 
state of matter through relativistic heavy ion collisions. This has led to experiments 
at BNL AGS and CERN SPS and to the building of the BNL Relativistic Heavy 
Ion Collider and a planning of the ALICE experiment at the CERN Large Hadron 
Collider. With the reported confirmations of the quark-hadron phase transition at 
the relativistic heavy ion collision experiments at the CERN SPS [[L|], the first step in 
the search for quark-gluon plasma, which pervaded the early universe, microseconds 
after the big bang and which may be present in the core of neutron stars, is complete. 

The emphasis of the experiments at the BNL RHIC and the CERN LHC will now 
necessarily shift to an accurate determination of the properties of the quark matter. 
An important observable for this is the speed of sound in the plasma, defined through: 

Often one writes, 

P = c 2 s e (2) 

for the equation of state of the quark matter, where e is the energy density and p is 
the pressure. For the simplest bag-model equation of state (with \i B = 0), we write 

e = 3aT 4 + B, (3) 
p = aT A — B, (4) 

7T 
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with B as the bag pressure and Nf as the number of flavours, so that c 2 s = 1/3. In 
general A = e — 3p measures the deviation of the equation of state from the ideal gas 
of massless quarks and gluons (when it is identically zero) and depends sensitively 
on the interactions present in the plasma. The lattice QCD calculations show that 
A > 0, till the temperature is several times the critical temperature JJ. This implies 
that in general c| < 1/3. Any experimental information on this will be most welcome. 

We show in the present work that the transverse momentum dependence of the 
survival probability of the J/ijj and T at RHIC and LHC energies are quite sensitive 
to the value of the speed of sound. The very long life-time of the plasma likely to be 
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attained at LHC makes it even more sensitive to the details of the equation of state 
of the quark matter through the transverse expansion of the plasma. 



2 Formulation 



The theory of quarkonium suppression in QGP || is very well studied and several 
excellent reviews exist ML [|, which dwell both on the phenomenology as well as 
on the experimental situation. We recall the basic details which are relevant for the 
present demonstration. 

The interquark potential for (non-relativistic) quarkonium states at zero temper- 
ature may be written as: 

V(r,0)=ar-- (6) 
r 

where r is the separation between Q and Q. The bound-states of cc and bb are 
well described if the parameters o = 0.192 GeV 2 , a =0.471, m c = 1.32 GeV, and 
rrib = 4.746 GeV are used 0. At finite temperatures the potential is modified due 
to colour screening, and evolves to: 



A*CT) 



a 



r 



The screening mass increases with temperature. When /x(T) — > 0, the equation (^) is 
recovered. At finite temperature, when r — > the 1/r behaviour is dominant, while 
as r —>■ oo the range of the potential decreases with fi(T). This makes the binding 
less effective at finite temperature. Semiclassically, one can write for the energy of 
the pair, 

-^( r > T) = 2rriQ H °—— + V(r, T) (8) 

rriQr z 

where < p 2 >< r 2 >= c = 0(oo). Radius of the bound state at any temperature 
is obtained by minimizing E(r,T). Beyond some critical value \xb for the screening 
mass /i(T), no minimum is found. The screening is now strong enough to make 
the binding impossible and the resonance can not form in the plasma. The ground 
state properties of some of the quarkonia reported by authors of Ref. are given in 
table 1. We have also listed the formation time of these resonances defined in Ref. || 
as the time taken by the heavy quark to traverse a distance equal to the radius of the 
quarkonium in its rest frame ~ m Q r QQ/PQQ-> where Pqq is the momentum of either of 
the quarks of the resonance. It may be recalled that somewhat different values for the 
formation time are reported by Blaizot and Ollitrault M who solve the bound-state 
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problem within the WKB approximation and define the formation time as the time 
spent by a quark in going between the two classical turning points. 

Now let us consider a central collision in a nucleus-nucleus collision, which results 
in the formation of quark gluon plasma at some time To- Let us concentrate at z = 
and on the region of energy density, e > e s which encloses the plasma which is dense 
enough to cause the melting of a particular state of quarkonium. We assume the 
plasma to cool, according to Bjorken's boost invariant (longitudinal) hydrodynamics 
and then generalize our results to include the transverse expansion of the plasma. 
We assume that the QQ pair is produced at the transverse position r at r = on 
the z = plane with momentum p T . In the collision frame, the pair would take a 



time equal to TpE^jM for the quarkonium to form, where = y Pt + M 2 and M 
is the mass of the quarkonium. During this time, the pair would have moved to the 
location (t + tfPt/M). If at this instant, the plasma has cooled to an energy density 
less than e s , the pair would escape and quarkonium would be formed. If however, the 
energy density is still larger than e s , the resonance will not form and we shall have a 
quarkonium suppression || [H| |TT] . 

It is easy to see that the px dependence of the survival probability will depend 
on how rapidly the plasma cools. If the initial energy density is sufficiently high, the 
plasma will take longer to cool and only the pairs with very high p? will escape. If 
however the plasma cools rapidly, then even pairs with moderate px will escape. The 
transverse expansion of the plasma can further accelerate the rate of cooling giving us 
an additional handle to explore the equation of state, which as we know, will control 
the expansion of the plasma. 



2.1 Longitudinal expansion of the plasma 

As indicated, we first take the Bjorken's boost-invariant longitudinal hydrodynamics 
to explore the expansion of the plasma. Thus, the energy momentum tensor of the 



plasma is written as [12 



= (e + p)u»u u + g^p, (9) 

where e is the energy-density, p is the pressure, and w M is the four velocity of the fluid, 
in a standard notation. If the effects of viscosity are neglected, the energy-momentum 
conservation is given by 

d^ v = 0. (10) 

The assumption of the boost-invariance provides that the energy density, pressure, 
and temperature become functions of only the proper time r and that the Eq.(|T0"D 
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simplifies to 



de 
dr 



e + p 



(11) 



The effect of the speed of sound is seen immediately. Using the Eq.(|2|), we can 
now write 



err 



l+d 



e(r )r 1+Cs = const. (12) 

so that if Cg is small, the cooling is slower. Chu and Matsui explored the con- 
sequence of the extremes cj. = 1/3 and c 2 s = on the pt dependence of the survival 
probability. We shall explore the sensitivity of the quarkonium suppression to the 
equation of state by (some-what arbitrarily) choosing two values of the speed of 
sound, l/y/3 and l/y/E, in the following. 

We now have all the ingredients to write down the survival probability and we 
closely follow Chu and Matsui for this. 

We take a simple parametrization for the energy-density profile: 



e(r ,r) 



r 

R 2 



6{R-r) 



(13) 



where r is the transverse co-ordinate and R is the radius of the nucleus. One can 
define an average energy density < eo > as 



nR 2 < e >= J 2nrdre(r) 



(14) 



so that 



e = (l + (3)<e > . (15) 

We have taken (3 = 1/2, which may be thought as indicative of the energy deposited 
being proportional to the number of participants in the system. In case one feels that 
the energy deposited may be proportional to number of nucleon-nucleon collisions 
then one can repeat the calculations with (3 = 1, which will reflect the proportional- 
ity of the deposited energy to the nuclear thickness. The average energy-density is 
obtained from the Bjorken formula: 

1 dE T 



< e >= 



7iRj>tq dy 



(16) 



where Et is the transverse energy deposited in the collision. 

The time r s when the energy density drops to e s is easily estimated as 



r s (r) 



To 



e(Vo,r) 
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T 
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/3/(l+c?) 
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As discussed earlier [10], we can equate the duration of screening r s (r) to the 



formation time tp = ^Tp for the quarkonium to get the critical radius, r s : 



R 



- / 7TF x (i+c*)//r 


1/2 
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V r s0 / 




r s0 - 



(18) 



where t s q = r s (r = 0). This critical radius, is seen to mark the boundary of the region 
where the quarkonium formation is suppressed. As discussed earlier, the quark-pair 
will escape the screening region (and form quarkonium) if its position and transverse 
momentum p T are such that 

|r + t fPt /M| > r s . (19) 

Thus, if is the angle between the vectors r and P t, then 

cos0 > [{r 2 s -r 2 )M - t 2 fP 2 t /m\ / [2r t fPt } , (20) 

which leads to a range of values of <fi when the quarkonium would escape. We also 
realize that if the right hand side of the above equation is greater than 1, then no 
angle is possible when the quarkonium can escape. Now we can write for the survival 
probability of the quarkonium: 

i-R 



'T 



R r+<t>max 
rdr 

11 "'- < ? i max 



.P(r,p T ) 



/ 



f 

2tt r dr P(r, P t) 
J o 



(21) 



where 0max is the maximum positive angle (0 < < tc) allowed by Eq.fl20D, and 

if 2/ < -1 



^max 



7T 

COS" 





\y\ 



(22) 



if -1< y < 1 , 
if y > 1 

where 

y=[( r 2 -r 2 )M-T 2 p 2 T /M] /[2rr F p T ] , (23) 

and P is the probability for the quark-pair production at r with transverse momentum 
p T , in a hard collision. Assuming that the P t and r dependence for hard collisions 
factor out, we approximate 

P(r, pr) = P(r, p T ) = f (r)g(p T ), (24) 

where we take 



f(r) oc 



1-1- 

R 2 



6{R-r) 



(25) 



with a = 1/2. The Eq.(pT|) can be solved analytically for some limiting cases of pp 
etc., see Ref. |IT| . 
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3 Transverse expansion of the plasma 



It is generally accepted that the rare-faction wave from the surface of the plasma will 
reach the centre by r = R/c s . For the case of lead nuclei, this comes to about 12 
fm/c. If the life time of the QGP is comparable to this time, the transverse expansion 
of the plasma can not be ignored. The transverse expansion of the plasma will lead 
to a much more rapid cooling than suggested by a purely longitudinal expansion. 



For the four-velocity of the collective flow we write: 

m m = (7»7v) 



(26) 



where v is the collective flow velocity and 7 = — v 2 . We further assume that 

the longitudinal flow of the plasma has a scaling solution, so that the boost-invariance 
along the longitudinal direction remains valid. Assuming cylindrical symmetry, valid 
for central collisions, it can be shown that the four-velocity should have the form, 



= 7 r (r, r)(t/r, v r cos0, v r sin0, z/t), 



(27) 



with 



T 

V 



1 -v 2 (r,r) 

(t 2 -z 2 ) 1 ' 2 

1. t + z 
- In . 

2 t-z 



-1/2 



(28) 



Thus all the Lorentz scalars are now functions of r and r, and independent of the 
space-time rapidity rj. This reduces the (3 + 1) dimensional expansion with cylindrical 
symmetry and boost-invariance along the longitudinal direction to: 



and 

where 

and 



d T T 00 + r^d r {rT m ) + r^T 00 + p) = 

d T T 01 + r- l d r [r(T 00 + p)v 2 ] + r~ l T m + d r p = 
T 00 = (e + p)u°u° - p 



T 01 = (e + p)u u u\ 



o„,i 



(29) 

(30) 
(31) 
(32) 



We see that the speed of sound which will appear through the dependence of the 
pressure p on the energy-density will affect the time-evolution of all the quantities, 
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especially through the gradient terms. This is of-course extensively documented. We 
solve these equations using well established methods JES] with initial energy density 
profiles as before and estimate the constant energy density contours [[TJ], [T5j appro- 
priate for e = e s to get r s (r). Rest of the treatment follows as before. In these 
calculations we have assumed the initial transverse velocity to be identically zero. 

We only need to identify the initial conditions. We consider Pb + Pb collisions 
(Au + Au, for RHIC) with the initial average energy densities: 



' 6.3 GeV/fm 3 SPS, r = 0.5 fm 
< eo >= J 60 GeV/fm 3 RHIC, r = 0.25 fm 



425 GeV/fm 3 LHC, r = 0.25 fm 



(33) 



The estimate for SPS is obtained from assumption of QGP formation in Pb + Pb 
experiments, while those for RHIC and LHC are taken from the self-screened parton 



cascade calculation |Uj. If a formation time for the SPS energies is assumed to be of 
the order of 1 fm/c, the estimate given above will drop to about 3 GeV/fm 3 . A larger 
initial energy density than the one assumed here for RHIC and LHC could be obtained 
by using the concept of parton saturation [17]. We are, however, interested in only a 
demonstration of the effect of equation of state on the quarkonium suppression, and 
we feel that the values used here are enough for this. 



4 Results 

4.1 Speed of sound vs. transverse expansion 

It is quite clear that a competition between the speed of sound and the onset of 
the transverse expansion during the life-time of the deconfining matter can lead to 
interesting possibilities. In order to illustrate this and to explore the consequences, we 
show in Fig. 1 the time corresponding to the constant energy density contours which 
enclose the deconfining matter- which can dissociate the directly produced J/ip (see 
table 1), at RHIC energies. We see that owing to the (relatively) short time that the 
QGP would take to cool down to ef'^, the effect of the transverse flow is marginal 
and for c 2 s = 1/3, limited to large radii. Large changes in the contour are seen when 
the speed of sound is varied. This is very important indeed. Note that the cooling 
to the value appropriate for T suppression is attained too quickly to be affected by 
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the transverse expansion, and even the change due to variation in the speed of sound 
is quite small. Of-course the duration of the deconfining medium is prolonged if the 
speed of sound is reduced. 

The corresponding results for the LHC energies are shown in Fig. 2. Now we see 
that at r = the duration of the deconfining medium reduces by a factor of 2 when 
the speed of sound is 1/ v 7 ^, and the transverse expansion of the plasma is allowed. 
This is a consequence of the longer time which the plasma takes to cool at LHC. 

The scenario for the T dissociating matter at LHC is quite akin to the case of J / ip 
dissociating medium at RHIC; the results are affected by the speed of sound and not 
by the transverse flow (Fig. 3). 

4.2 Consequences for survival of quarkonia 

Now we return to the transverse momentum dependence of the survival probabilities. 
As a first step, we plot the survival of the directly produced J/^'s at SPS, RHIC 
and LHC energies (Fig. 4) when only longitudinal expansion is accounted for and 
the speed of sound is varied. We see that RHIC energies provide the most suitable 
environment to measure the speed of sound with the help of J/ip suppression. The 
variations in the pt dependence is too meagre at SPS (due to a very short duration 
of the deconfining medium) and at LHC (now, due to a very long duration!) when 
the speed of sound is varied. 

This advantage of RHIC energies is maintained when the transverse expansion is 
accounted for (Fig. 5). As one could have expected from the contours (Fig. 1), the 
results are more sensitive to the variation of the speed of sound than to the transverse 
flow. The accuracy of the procedure is seen from the fact that the survival probability 
around pt = 15 GeV for c 2 s = 1/3 is identical for the longitudinal and the transverse 
expansion of the plasma. This is a direct reflection of the identity of the corresponding 
contours near r = (Fig. 1). 

The J/ip suppression at LHC energies, as indicated, becomes sensitive to the 
transverse flow, the shape of the survival probability changes and the largest pr for 
which the formation is definitely possible is enhanced by about 10 GeV (Fig. 6). 

The T suppression, which in our prescription is possible only at the LHC energy, is 
seen to be clearly affected by the speed of sound but not by the transverse expansion 
of the plasma (Fig. 7). 



9 



4.3 Consequences of chain decays of quarkonia 



The entire discussion so far has been in terms of the directly produced J/ip's and 
T's. However it is well established that only about 58% of the observed J /if) in pp 
collisions originate directly, while 30% of them come from \ c decay and 12% from the 
decay of ip'. Thus the survival probability of the J/ip in the QGP can be written as: 

S = 0.585^ + 0.3 S Xc + 0.125^ , (34) 

in an obvious notation. We give the survival probabilities of these resonances for a 
transversely expanding plasma at RHIC with the speed of sound as 1/ \/3 in the left 
panel of Fig. 8. We see that the competition of the formation times and the duration 
of the sufficient dissociation energies render a rich detail to the suppression pattern of 
the charmonia. We shall later argue that the ifi n s are easily destroyed by a moderately 
hot hadronic gas as well, as their binding energy is on the order of just 50 MeV. The 
right panel of the figure shows the survival probability as a function of the transverse 
momentum when the speed of sound in the plasma is varied from l/v3 to 1/y/E. We 
see that while the gross features for shape of the survival probability remain similar 
to Fig. 5, as seen earlier, the survival of the J/ip for larger p? is now enhanced as the 
Xc's, which decay to form J/ip start escaping. Thus the complete escape of the Xc 
having px > 10 GeV (for cj. = 1/3) lends a distint kink in S(pr) at Pt ~ 10 GeV. 
Its location, which shifts to about 14 GeV when the speed of sound is decreased, can 
perhaps be more accurately determined by plotting the derivative dS/dpx, which will 
have a discontinuity there. If the statistics is really good (which unfortunately is a 
somewhat difficult proposition), this discontinuity in the dS/dpT for the J/ip can be 
a unique signature of melting of the resonance in plasma, as the pt dependence of the 
survival probability due to absorption by hadrons should be weak |0| and smooth ||. 



The corresponding results for the LHC energies are given in Fig. 9, in an analogus 
manner, and we see the characteristic 'kink' in S(px) has shifted to about p? = 16 
GeV, announcing the complete escape of Xc from the plasma. 

The decay contribution of the resonances completely alters the shape of the sur- 
vival probabilities for T from that seen earlier (Fig. 7). 

We note that (Fig. 10) both Xb and T' get suppressed at the RHIC energies, 
while the directly produced T is likely to escape. However, as only about 54% of 
the T's may be directly produced, while about 32% have their origin in the decay of 
Xb and a further 14% in the decay of T', the resultant (right panel, Fig. 10) survival 
probability is just about 50% for the lowest p?, signalling the suppressions of the 
higher resonances. 
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The results for LHC energies become quite dramatic (Fig. 11) as now the 'kink' 
in the survival probability is very clearly seen at px ~20 GeV for c 2 s = 1/3 and at 
~ 26 GeV for the lower speed of sound. If one looks at the dS/dpx then it could be 
very useful indeed. We may add that fluctuations in the initial conditions etc. may 
perhaps make it difficult to notice this aspect for J/ip's, though for the T's these 
should survive, provided we have good statistics at these px- 

The large difference in this behaviour seen between charmonium and bottomonium 
suppression has its origin in the large difference in the energy densities required to 
melt the different resonances of charmonia and bottomonia. 



4.4 Absorption by nucleons and comovers 

So far we have discussed the fate of quarkonia only when the presence of quark 
gluon plasma is considered. It is very well established that there are several aspects 
like initial state scattering of the partons, shadowing of partons, absorption of the 
pre-resonances ( \QQg > states) by the nucleons before they evolve into physical 
quarkonia, and also dissociation of the resonances by the comoving hadrons |6| . It 
has been argued that the absorption by co-moving hadrons will be important for if)', 
due to its very small binding energy, while for more tightly bound resonances it may 
be weak [Q. 

Let us briefly comment on them one-by-one. Shadowing of partons should play 
an important role in the reduced production of quarkonia, especially at the LHC 
energies. It is clear that if shadowing is important, we shall witness a larger effect on 
J /if} than on T, because of the smaller values of the x for gluons. At the same time, 
the effect of shadowing should be similar for different resonances of the charmonium 
(or bottomonium), as similar x values would be involved for them. 

The absorption of the pre-resonances by the nucleons is another source of pr 
dependence. It is important, to recall once again that as the absorption is operating 
on the pre-resonance, the effect should be identical for all the states of the quarkonium 
which are formed. 

This is a very important consideration as it is clear that if we look at the ratio 
of rates for different states of J/if> or the T family as a function of pt, then in the 
absence of QGP-effects they would be identical to what one would have expected 
in absence of nuclear absorption and shadowing, providing a clear pedestal for the 
observation of QGP |TS] . 
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There is another aspect of pt dependence which needs to be commented upon, 
before we conclude. The (initial state) scattering of partons, before the gluons of the 
projectile and the target nucleons fuse to produce the QQ-pair, leads to an increase of 
the < p T > of the resonance which emerges from the collision The increase in the 
< p\ >, compared to that for pp collisions is directly related to number of collisions 
the nucleons are likely to undergo, before the gluonic fusion takes place. This leads 
to a rich possibility of relating the average transverse momentum of the quarkonium 
to the transverse energy deposited in the collision (which decides the number of 
participants and hence the number of collisions). Considering that collisions with 
large Et may have formation of QGP in the dense part of the overlapping region, the 
quarkonia, which are produced in the densest part (and hence contributing the largest 
increase in the transverse momentum) are also most likely to melt and disappear. 
This may lead to a characteristic saturation and even turn-over of the < p T > when 
plotted against Et when the QGP formation takes place. In absence of QGP, this 
curve would continue to rise with Et- 

Obviously, all these (well explored and yet non-QGP) effects need to be accounted 
for, before we can begin to see the suppression of the quarkonium due to the formation 
of QGP. It seems that this has been achieved at least at the SPS energies P, §]. 

The next step would obviously be the one discussed in the present work, that of 
looking for the pt dependence of the survival probability, to see if we can get more 
detailed information on the equation of the state of the plasma. This would require 
high precision data, extending to larger pt at several Et- This may prove to be 
difficult, though not impossible in principle at least. It will however prove to be very 
valuable, if it can be done. 

5 Summary and Discussion 

We have seen that the survival probability of J/ip at RHIC energies and that of T 
at LHC energies can provide valuable information about the equation of state of the 
quark matter, as the results are not affected by uncertainties of transverse expansion 
of the plasma. If the transverse expansion of the plasma takes place, it gives a distinct 
shape to the survival probability of the J/ip at LHC energies, whose detection will 
be a sure signature of the transverse flow of the plasma within the QGP phase. 

Before concluding we would add that in these exploratory demonstrations we 
have chosen some specific values [§] for the deconfining matter which can dissociate 
quarkonia. These could be different, in particular the e J J^ could be much larger than 
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the value used here. This, however, will not change the basic results as the time-scales 
involved in the T and the J/ijj suppression are so very different. 

The other uncertainty comes from the usage of e s as the criterion for deconfine- 
ment. One could have as well used the Debye mass as the temperature and the 
fugacity changed in a chemically equilibrating plasma, to fix the deconfining zone. 
This can indeed be done, along with the other extreme of the Debye mass estimated 
from lattice QCD. This has been studied in great detail by authors of Ref. |H| and 
can be easily extended to the present case. We plan to do it in a future publication. 

In brief, we have shown that the J/ip and T suppression at RHIC and LHC 
energies can be successfully used to map the equation of state for the quark-matter. 
As the two processes will map different but over-lapping regions, taken together, these 
results will help us to explore a vast region of the equation of state. 

We thank Bikash Sinha for useful discussions. 
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Table 1: Critical screening masses, etc. for quarkonia |7|, [^]. 



State 




Xc 


if/ 


T 




Xb 


M (GeV) 


3.1 


3.5 


3.7 


9.4 


10.0 


9.9 


r (fm) 


0.45 


0.70 


0.88 


0.23 


0.51 


0.41 


r F (fm) 


0.89 


2.0 


1.5 


0.76 


1.9 


2.6 


Hd (GeV) 


0.70 


0.34 


0.36 


1.57 


0.67 


0.56 


T d /T c 


1.17 


1.0 


1.0 


2.62 


1.12 


1.0 


e s (GeV/fm 3 ) 


1.92 


1.12 


1.12 


43.37 


1.65 


1.12 
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Life-time of J/Y'-dissociating deconfined matter 




2 4 6 

r (fm) 

Figure 1: The contour of deconflning matter capable of suppressing direct J/ip for- 
mation at RHIC energies. The solid curves are results with transverse flow, while the 
dashed curves are for longitudinal flow. 
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Figure 2: The contour of deconfming matter capable of suppressing direct J/ip for- 
mation at LHC energies. The solid curves are results with transverse flow, while the 
dashed curves are for longitudinal flow. 
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Life-time of T-dissociating deconfined matter 
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Figure 3: The contour of deconfining matter capable of suppressing direct T formation 
at LHC energies. The solid curves are results with transverse flow, while the dashed 
curves are for longitudinal flow. 
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Survival of directly produced J/^ 
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Figure 4: The survival probability for J/ip at SPS, RHIC, and LHC energies, with 
only longitudinal cooling, when speed of sound is changed. 
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Survival of directly produced J/ip 
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Figure 5: The survival probability for J/ip at RHIC energy with longitudinal and 
transverse cooling when speed of sound is changed. 
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Survival of directly produced J/ij/ 
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Figure 6: The survival probability for J/ip at LHC energy with longitudinal and 
transverse cooling when speed of sound is changed. 
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Survival of directly produced T 
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Figure 7: The survival probability for T at (RHIC and) LHC energy with longitudinal 
and transverse cooling when speed of sound is changed. 
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Charmonia, Au+Au@RHIC J/^, direct+decay 

i i i i i i i i i i i i i i i i i i i i i i i i 1 1 1 1 i i i 1 1 i i i i i 1 1 1 1 1 i 1 1 1 1 1 




p T (GeV) p T (GeV) 

Figure 8: The survival probability for directly produced J ftp, Xc and ip' at RHIC 
energies, for a transversely expanding plasma (left panel), for c 2 s = 1/3. The right 
panel shows the results when the decays of the resonances is accounted for. 
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Charmonia, Pb+Pb@LHC 



direct+decay 




Bottomonia, Au+Au@RHIC 



T, direct+decay 




Figure 10: The survival probability for directly produced T, \c and T' at RHIC 
energies, for a transveresly expanding plasma (left panel), for c 2 s = 1/3. The right 
panel shows the results when the decays of the resonances is accounted for. 
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Bottomonia, Pb+Pb@LHC 



T, direct+decay 




